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Abstract 

We apply the perturbative chiral quark model (PCQM) at one loop to analyse 
the electromagnetic structure of nucleons. This model is based on an effective 
Lagrangian, where baryons are described by relativistic valence quarks and a 
perturbative cloud of Goldstone bosons. Including the electromagnetic inter- 
action we first develop the formalism up to one-loop in the Goldstone boson 
fluctuation relying on renormalization by use of counterterms. Local gauge 
invariance is satisfied both on the Lagrangian level and also for the relevant 
baryon matrix elements in the Breit frame. We apply the formalism to ob- 
tain analytical expressions for the nucleon charge and magnetic form factors, 
which are expressed in terms of fundamental parameters of low-energy pion- 
nucleon physics (weak pion decay constant, axial nucleon coupling, strong 
pion-nucleon form factor) and of only one model parameter (radius of the 
nucleonic three-quark core). A detailed numerical analysis for the nucleon 
magnetic moments, charge and magnetic radii and also for the momentum 
dependence of form factors is presented. 
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I. INTRODUCTION 

The study of the electromagnetic form factors of the nucleon is of fundamental im- 
portance in hadron physics to gain a deeper understanding of the baryon structure and 
the interplay between strong and electromagnetic interactions. Completed and ongoing 
experiments at ELSA, JLab, MAMI, MIT, NIKHEF and other laboratories on improved 
measurements of the elastic proton and neutron form factors stimulated a comprehensive 
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theoretical study of these quantities. The theoretical description of electromagnetic form 
factors was performed in detail within different approaches of low-energy hadron physics: 
QCD Sum Rules, Chiral Perturbation Theory, relativistic and non-relativistic quark mod- 
els, QCD-motivated approaches based on solutions of the Schwinger-Dyson/Bethe-Salpeter 
and Faddeev equations, soliton-type models, etc. For recent experimental and theoretical 
advances with respect to the electromagnetic structure of nucleons see the proceedings of 
the last conferences |[J- H. 

Since the early eighties chiral quark models ||, describing the nucleon as a bound 
system of valence quarks with a surrounding pion cloud, play an important role in the 
description of low-energy nucleon physics. These models include the two main features 
of low-energy hadron structure, confinement and chiral symmetry. With respect to the 
treatment of the pion cloud these approaches fall essentially into two categories. 

The first type of chiral quark models assumes that the valence quark content dominates 
the nucleon, thereby treating pion contributions perturbatively • Originally, this idea 

was formulated in the context of the cloudy bag model ||. By imposing chiral symmetry 
the MIT bag model [[J was extended to include the interaction of the confined quarks with 
the pion fields on the bag surface. With the pion cloud treated as a perturbation on the 
basic features of the MIT bag, pionic effects generally improve the description of nucleon 
observables. Later, similar perturbative chiral models were developed where the 

rather unphysical sharp bag boundary is replaced by a finite surface thickness of the quark 
core. By introducing a static quark potential of general form, these quark models contain a 
set of free parameters characterizing the confinement (coupling strength) and/or the quark 
masses. The perturbative technique allows a fully quantized treatment of the pion field up 
to a given order in accuracy. Although formulated on the quark level, where confinement 
is put in phenomenologically, perturbative chiral quark models are formally close to chiral 
perturbation theory on the hadron level. 

Alternatively, when the pion cloud is assumed to dominate the nucleon structure this 
effect has to be treated non-perturbatively. The non-perturbative approaches are based for 
example on Refs. 0, where the chiral quark soliton model was derived. This model is based 
on the concept that the QCD instanton vacuum is responsible for the spontaneous breaking 
of chiral symmetry, which in turn leads to an effective chiral Lagrangian at low energy as 
"derived" from QCD. The classical pion field (the soliton) is described by a trial profile 
function, which is fixed by minimizing the energy of the nucleon. Further quantization of 
slow rotations of this soliton field leads to a nucleon state, which is built up from rotational 
excitations of the classical nucleon. On the phenomenological level the chiral quark soliton 
model tends to be advantageous in the description of the nucleon spin structure, that is 
for large momentum transfers, but is comparable to the original perturbative chiral quark 
models in the description of low-energy nucleon properties. 

As a further development of chiral quark models with a perturbative treatment of the 
pion cloud @-|6|,|8[], we extended the relativistic quark model suggested in || for the study 
of the low-energy properties of the nucleon ]nj . Compared to the previous similar models 
of Refs. our current approach contains several new features: i) generalization of the 

phenomenological confining potential; ii) SU(3) extension of chiral symmetry to include the 
kaon and eta-meson cloud contributions; iii) consistent formulation of perturbation theory 
both on the quark and baryon level by use of renormalization techniques and by allowing 
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to account for excited quark states in the meson loop diagrams; iv) fulfillment of the con- 
straints imposed by chiral symmetry (low-energy theorems), including the current quark 
mass expansion of the matrix elements (for details see Ref. ||10|| ); v) possible consistency 
with chiral perturbation theory as for example demonstrated |H| for the chiral expansion 
of the nucleon mass. In the following we refer to our model as the perturbative chiral quark 
model (PCQM). 

The PCQM is based on an effective chiral Lagrangian describing quarks as relativistic 
fermions moving in a self-consistent field (static potential). The latter is described by a 
scalar potential S providing confinement of quarks and the time component of a vector 
potential responsible for short-range fluctuations of the gluon field configurations [11 



(see also recent lattice calculations |12|). The model potential defines unperturbed wave 
functions of quarks which are subsequently used in the calculation of baryon properties. 
Baryons in the PCQM are described as bound states of valence quarks surrounded by a 
cloud of Goldstone bosons (n,K,ri) as required by chiral symmetry. Interaction of quarks 
with Goldstone bosons is introduced on the basis of the nonlinear a- model [ 13| . When 



considering mesons fields as small fluctuations we restrict ourselves to the linear form of the 
meson-quark interaction. With the derived interaction Lagrangian we do our perturbation 
theory in the expansion parameter 1/F (where F is the pion leptonic decay constant in the 
chiral limit). We also treat the mass term of the current quarks as a small perturbation. 
Dressing the baryon three-quark core by a cloud of Goldstone mesons corresponds to the 
inclusion of the sea-quark contribution. All calculations are performed at one loop or at 
order of accuracy o(l/F 2 ,rh,m s ). The chiral limit with m,m s — > is well defined. 

In Ref. |T(J our approach was successfully applied to the study of sigma-term physics. Our 
result for the ttN sigma term u^n ~ 45 MeV is in good agreement with the value deduced 
by Gasser, Leutwyler and Sainio ]TJ| using dispersion-relation techniques and exploiting the 
chiral symmetry constraints. A complete analysis of meson-nucleon and meson-delta-isobar 
sigma-terms, is given including the strangeness content of the nucleon. We have a physically 



reasonable prediction for the quantity = 0.076 ± 0.012 |10[ characterizing the strange 
quark contribution to the nucleon mass. Our predictions for the strangeness content of the 
nucleon and for KN sigma-terms are important for the ongoing DA$NE experiments at 



Frascati |L5] 



In this paper, we apply the PCQM to analyse the electromagnetic structure of nucleons. 
The main intention of the present work is to develop the full formalism for the effective 
Lagrangian of the PCQM, where Goldstone boson fluctuations are treated as higher order 
effects to the pure valence quark core results. We give detailed consideration of gauge 
invariance in noncovariant chiral quark model treating pion cloud perturbatively. 

In the present article we proceed as follows. In the following section, we first describe the 
basic notions of our approach: the underlying effective Lagrangian, choice of parameters and 
fulfilment of low-energy theorems. We extend our model Lagangian by including the electro- 
magnetic interaction. A consistent formalism of the perturbation series up to one meson loop 
and up to terms linear in the current quark masses is presented. We utilize a renormalization 
technique, which, by introducing counterterms, greatly simplifies the evaluation. Thereby, 
local gauge invariance of the electromagnetic interaction is shown to be fulfilled both on 
the Lagrangian level and for baryon matrix elements set up in the Breit frame. In Sect. 3 
we concentrate on the specific application of the PCQM to the electromagnetic properties 
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of the nucleoli. We derive analytical expressions for the nucleon charge and magnetic form 
factors expressed in terms of fundamental parameters of low-energy pion-nucleon physics 
(weak pion decay constant, axial nucleon coupling, strong pion-nucleon form factor) and 
of only one model parameter (radius of the nucleonic three-quark core). Numerical results 
in comparison with data are presented to test the basic and standard phenomenological 
implications of the model. Finally, Sect. 4 contains a summary of our major conclusions. 



II. PERTURB ATIVE CHIRAL QUARK MODEL 



A. Effective Lagrangian and zeroth order properties 



The perturbative chiral quark model (PCQM) ]TUj is based on an effective chiral La- 
grangian describing the valence quarks of baryons as relativistic fermions moving in a self- 
consistent field (static potential) V e ff(r) = S(r) + r y°V(r) with r = |x| which are 
supplemented by a cloud of Goldstone bosons (n,K,ri). Treating Goldstone fields as small 
fluctuations around the three-quark (3q) core we derive a linearized effective Lagrangian 
C e ff. The Lagrangian C e ff 
invariant term 



Cy™ v + £-xSB, derived in Ref. [I(J, includes the linear chiral- 



= - S(r) - 7 t V(r)]V>(z) + -(fy$(a:)) 2 - ^(x)S(r)ir 



i(}(x) 



and a term C x sb which explicitly breaks chiral symmetry 

C x sb(x) = -i>(x)Mi/>(x) - f Tr{$(x), {$(x), M}} 



(2) 



containing the mass terms for quarks and mesons |10| . The octet matrix $ of pseudoscalar 
mesons is defined as: 



E 



V2 U V2 



\ 



71 

K~ 



-7T 



7T+ K+ \ 

°/V2 + v /Vq k° 

K° -2rj/y/6 ) 



(3) 



F = 88 MeV [p~Q] , p~T[| is the pion decay constant in the chiral limit; A4 = diag{m, m, m s } 
is the mass matrix of current quarksQ B = — < 0\uu\0 > / F 2 is the low-energy constant 
which measures the vacuum expectation value of the scalar quark densities in the chiral 



limit 1 18]. We rely on the standard picture of chiral symmetry breaking and for the 
masses of pseudoscalar mesons we use the leading term in their chiral expansion (i.e. linear 
in the current quark mass): 



Ml = 2mB, 



M 2 K = (rh + m s )B, 



Ml 



irk + 2m s )B. 



(4) 



Here we restrict to the isospin symmetry limit with m u = = rh. 



4 



Meson masses obviously satisfy the Gell-Mann-Oakes-Renner (|j) and the Gell-Mann-Okubo 

4M](. In the evaluation we use the following set of QCD parame- 



relation 3M 2 + Ml 



ters 19 



m 



7 MeV, 



m 



25 and B 



Ml 



2rh 



1.4 GeV. 



(5) 



consis- 



Furthermore, the linearized effective Lagrangian fulfils the PCAC requirement |L0 
tent with the Goldberger-Treiman relation (as shown later on). 

To derive the properties of baryons, which are modelled as bound states of valence 
quarks surrounded by a meson cloud, we formulate perturbation theory. At zeroth order, 
the unperturbed Lagrangian simply describes the nucleon by three relativistic valence quarks 
which are confined by the effective one-body potential in the Dirac equation. The mass m c § re 
of the three-quark core of the nucleon is then related to the single quark ground state energy 
So by m c ^ re — 3 • S q (Here we do not discuss the removal of the spurious contribution to the 
baryon mass arising from the centre-of-mass motion of the bound state. It will be subject 
of a future publication.). For the unperturbed three-quark state we introduce the notation 
|0o > with the normalization < 0o|0o >= 1- The single quark ground state energy £ and 
wave function (w.f.) Uq(x) are obtained from the Dirac equation 



[-iaV + 0S(r) +V(r)- £ }u (x) 



0. 



(6) 



The quark w.f. u (x) belongs to the basis of potential eigenstates (including excited quark 
and antiquark solutions) used for expanding the quark field operator ip(x). Here we restrict 
the expansion to the ground state contribution with 



-0(x) = b u (x) exp(-i£ t), 



(7) 



where &o is the corresponding single quark annihilation operator. 

At the unperturbed level the current quark mass is not included in Eq. (|J) to simplify our 
calculational technique. Instead we consider the quark mass term as a small perturbation 
[T9] l, as will be discussed later on. Inclusion of a finite current quark mass leads to a 
displacement of the single quark energy which, for example, is relevant for the calculation 



of the sigma-term (see discussion in Ref. [l(J); a quantity which vanishes in the chiral limit. 
On the other hand, the effect of a finite current quark mass on observables which survive in 
the chiral limit, like magnetic moments, charge radii, etc., is quite negligible. 

For a given form of the potentials S(r) and V(r) the Dirac equation can be solved 



numerically. Here, for the sake of simplicity, we use a variational Gaussian ansatz |20| for 
the quark wave function given by the analytical form: 



u (x) = N exp 



x' 



2R 2 



where 



• ax I Xf Xc , 



N=[n 3 / 2 R 3 (l + 3p 2 /2)}- 1 / 2 
is a constant fixed by the normalization condition 



(8) 



(9) 
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d 3 xul(x) u (x) = 1; (10) 

Xs, X/; Xc ar e the spin, flavor and color quark wave functions, respectively. Our Gaussian 
ansatz contains two model parameters: the dimensional parameter R and the dimensionless 
parameter p. The parameter p can be related to the axial coupling constant qa calculated 
in zeroth-order (or 3q-core) approximation: 



9a= oU-T-nn^) = o^r— > ( n ) 



5/ 2p 2 \ _ 51 + 2 7 
3V 1 + fpV 3 

where 7 is a relativistic reduction factor 



1 - §p 2 9 1 

t = TT|7 = To^-2- < 12) 



Note that since PCAC is fulfilled in our model |L0| , on the tree level the axial charge g& is on 



the tree level related to the pion-nucleon coupling constant G^nn by the Goldberger-Treiman 
relation 

5m N ( 2p 2 \ _ m N 

- W I 1 - TTJ? ) = ~T 9A - (13) 

where is the physical nucleon mass. 

The parameter R can be physically understood as the mean radius of the three-quark 
core and is related to the charge radius < r\ >£ of the proton in the leading-order (or 
zeroth-order) approximation as [|l(J 

<r|> ^i±44.w 2 _iy (1 4, 



1 + fp 2 V 2, 

In our calculations we use the value £4=1.25 as obtained in ChPT []I7|. We therefore have 



only one free parameter, R. In the numerical evaluation R is varied in the region from 0.55 
fm to 0.65 fm corresponding to a change of < r\ >£ G in the region from 0.5 fm 2 to 0.7 
fm 2 . The use of the Gaussian ansatz ([8]) in its exact form restricts the scalar confinement 
potential S(r) to 

S( r ) = 1 ~ 3p2 + -E- r 2 = M + cr 2 , (15) 
V ; 2pR 2R^ ' v ; 

= M =c 

expressed in terms of the parameters R and p. The constant part of the scalar potential M 
can be interpreted as the constituent mass of the quark, which is simply a displacement of 
the current quark mass due to the potential S(r). The parameter c is a constant defining the 
radial (quadratic) dependence of the scalar potential. Numerically, for our set of parameters, 
we get M = 230 ± 20 MeV and c = 0.007 ± 0.002 GeV 3 . Error bars for theoretical results 
quoted here and in the following correspond to the variation of the parameter R. For the 
vector potential we get the following expression 

where the single quark energy Sq is a free parameter in the Gaussian ansatz. 
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B. Perturbation theory and nucleon mass 



Before setting out to present the renormalization scheme of the PCQM, we first define 
and discuss the quantities, relevant for mass and wave function renormalization. Following 
the Gell-Mann and Low theorem [21| we define the mass shift of the nucleonic three-quark 
ground state Amjv due to the interaction with Goldstone mesons as0 

oo -n „ 

Am N = ^< 0o| X — J iS(ti)d 4 x 1 ...d 4 x n T[£i(x 1 )...£ I (x n )}\ ( j) >? . (17) 

n=l 

In Eq. (|T7p the strong interaction Lagrangian Cj treated as a perturbation is defined asf] 

d(x) = -^(a;)i7 5 ^5(r)^(x) (18) 



and subscript "c" refers to contributions from connected graphs only. We evaluate Eq. (|Tj 
at one loop to order o(l/F 2 ) using Wick's theorem and the appropriate propagators. For the 
quark field we use a Feynman propagator for a fermion in a binding potential. By restricting 
the summation over intermediate quark states to the ground state we get 

iG^{x : y) =< (f) \T{i/j(x)i>(y)}\(j) >-> u (x)u (y) exp[-i£ (x - y o )}0(x o - y ). (19) 

For meson fields we use the free Feynman propagator for a boson field with 

i/Mx-v) =< 0|T{. |( ,).,M}|0 >= kj j^ J^-^-_f . (20) 



Superscript "iV" in Eq. (|17|) indicates that the matrix elements are projected on the re- 
spective nucleon states. The nucleon wave function |iV > is conventionally set up by the 
product of the SU(6) spin-flavor w.f. and SU(3) C color w.f. (see details in f23f), where the 
nonrelativistic single quark spin w.f. is replaced by the relativistic ground state solution of 
Eq. (H). Projection of "one-body" diagrams on the nucleon state refers to 



xWs> lffjS ' S XfXs ^<iV|X(/J)«|iV> (21) 

i=l 



where the single particle matrix element of the operators I and J, acting in flavor and spin 
space, is replaced by the one embedded in the nucleon state. For "two-body" diagrams with 
two independent quark indices i and j the projection prescription reads as 



2 We do not take into account the electromagnetic contribution to the nucleon mass since this 
effect is quite negligible. 

3 In the calculation of matrix elements we use the interaction Lagrangian and Wick's T-ordering 
product for field operators p^j . 
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X*f>xl>ll' f Jl' s XfXs ®xUl I^J^XkX* P -^<N\j:(I 1 J 1 f)®(I 2 J 2 )^\N>. (22) 

The total nucleoli mass is given by m r N = m c § re + Am^. Superscript "r" refers to the 
renormalized value of the nucleon mass at one loop, that is the order of accuracy we are 
working in. The diagrams that contribute to the nucleon mass shift Am^ at one loop are 
shown in Fig.l: meson cloud (Fig. la) and meson exchange diagrams (Fig. lb). The explicit 
expression for the nucleon mass including one-loop corrections is given by 



m r M = m core 



N - ,,o N + Am N = 3 {€ + 7 m) + £ d% n(M|) (23) 

<E>=7r,A",7/ 

wifh rT* — rl K — —rT* rf> - ^-rF 

Willi a N — , a N — Yg"^' a N ~ 57 a Ni 

where d% are the recoupling coefficients defining the partial contribution of the 7T, K and rj- 
meson cloud to the mass shift of the nucleon. For the following it is also useful to separately 
indicate the contributions to d% from the meson cloud (drf ) and the meson exchange 
diagrams (dj ): 

d N - 40Q , d N -d N - m , d N - 40Q (24 ) 



,tt;EX _ 90 jK;EX _ n J^EX _ 6 

J JV ~~ /inn' ^iV — U, Ctjy 



400' iV ' JV 400 

The self-energy operators I1(M|), corresponding to meson cloud contributions with definite 
flavor, differ only in their value for the meson mass and are given by 

For a meson with three- momentum p the meson energy is w<$,(p 2 ) = J M| + p 2 with p = \p\ 
and F-nNNip 2 ) is the ttNN form factor normalized to unity at zero recoil (p 2 = 0): 

W^-^.^-A)}. (26) 

Finally, the effect of a finite current quark mass rh on the nucleon mass shift is taken into 
account perturbatively (for details see [TO), resulting in the linear term 37m in Eq. (|23|). 



C. Renormalization of the PCQM 



To redefine our perturbation series up to a given order in terms of renormalized quantities 
a set of counterterms SC has to be introduced in the Lagrangian. Thereby, the counterterms 
play a dual role: i) to maintain the proper definition of physical parameters, such as nucleon 
mass and, in particular, the nucleon charge and ii) to effectively reduce the number of 
Feynman diagrams to be evaluated. 



S 



1. Renormalization of the quark field 



First, we introduce the renormalized quark field ip T with renormalized mass A4 r , substi- 
tuting the original field ip. Again, we restrict the expansion of the renormalized quark field 
to the ground state with 



W(x; m\) = b u r (x; m\) exp(-i£J(m[)t), 



(27) 



where "z" is the SU(3) flavor index; £o(m£) is the renormalized energy of the quark field in 
the ground state obtained from the solution of the Dirac equation 



-iaV + pml + pS(r) + V{r) - £ r (m r )K(f; m r ) = 0. 



(28) 



Using the derivations of the previous chapter, the renormalized mass m[ of the quark field 
is given by 



m, t 



m 



m — 6m 



(29) 



6m o 



m H + — 

3 7 



d 



^ MC U(M 2 K ) + 2dl MC 



The meson exchange contribution will be included when introducing nucleon mass renor- 
malization. For the quark masses we will use in the following the compact notation: 



A4 r = diag{m r , m r , m r s } and 6A4 = diag{<5m, 6m, Sm s }. 



(30) 



The solutions of Eq. (|28"D, Sq {m\) and u r (x; ml), are functions of m\. Obviously, the differ- 
ence between nonstrange and strange quark solutions is solely due to the flavor dependent 
quark mass m^. In the limit m\ — > the solutions for nonstrange and strange quarks are 
degenerate: E^rrf) = So and Uq(x; 0) = uo(x) (see Eq. (||)). For the renormalized wave 
function iiQ(x;ml) we again consider the Gaussian ansatz 



u {x;m it 



with normalization 



N(ml) exp 



-c(m[) 



x ■ 



2R 2 



ip(m r t ) ^ 



d A x Uq(x; ml) u r Q {x; m[) = 1. 



Xs Xf Xc 



(31) 



(32) 



In Eq. ( 
follows: 



the functions iV(m[), c(m[) and p(m[) are normalized at the point m\ = as 



N(0) = N, 



c(0) 



p(0) = p 



(33) 



The product p(m[)c(m[) can be shown to be m[-invariant and we therefore obtain the 
additional condition 



p(m[) c(m\) = p. 



(34) 
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Treating TTIa clS cL small perturbation, Eq. (B3) can be solved perturbatively, resulting in: 

^o( m i ) = £o + ^o( m D and «o(x; m [) = ^(a 7 ) + 5u (x; m[) (35) 

where 



dSoim'i) = 7m- and 5uQ(a;;m[ 



ml pi? ^ + TP 2 



2 1 + f p 2 V 1 + |p 2 i? 2 



-7^ + 7 «oW- 



For our set of model parameters the ground state quark energy £q is about 400 MeV and 
for the energy corrections 6£q relative to £q we obtain 



5£o(rh r 



£0 



14% and 



S£o(m r . 



£0 



/o. 



(36) 



Given the small corrections expressed in Eq. (p6|), the perturbative treatment of a finite 
(renormalized) quark mass is a meaningful procedure. 



2. Renormalized effective Lagrangian 

Having set up renormalized fields and masses for the quarks we are in the position to 
rewrite the original Lagrangian. The renormalized effective Lagrangian including the photon 
field is now written as 

£- r fuii = £ip + + £ph + £int- (37) 
The renormalized quark Lagrangian £^ defines free nucleon dynamics at one-loop with 

C; = % + Cfap, % = V{x)[i @-M r - S(r) - 7 V(r)]^(x) (38) 

i=l i=4 

The parameters 5A4 of Eq. (|30"D guarantee the proper nucleon mass renormalization due 
to the meson cloud diagrams of Fig. la. The terms contained in £/^ a are introduced for 
the purpose of nucleon mass renormalization due to the meson exchange diagram of Fig. lb. 
The corresponding renormalization parameters c n , ck and c v are deduced from Eqs. (p3|) 
and (121) as 



c *--4T^ n(M|) ' (39) 



The free meson Lagrangian C<$> is written as 



= ~\ E + M^(x) (40) 



where □ = d^d^ and M?- is the diagonal meson mass matrix with 
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Ml = M 2 22 = M 3 2 3 = Ml 



Ml = Ml = Ml, = M 7 2 7 = M|, 



Mi 8 = M>. (41) 



For the photon field we have the usual kinetic term 

C ph = ~F^(x)F^(x) with F lu/ (x)=d v A li (x)-d ll A v (x). (42) 



The renormalized interaction Lagrangian C int = C str + C r em contains a part due to the strong 

C r str = Cf + SC str (43) 

and the electromagnetic interaction 

C em = Cj m + 5C em . (44) 

The strong interaction term Cf is given by 



Cf = -^(x)t^^^S{r)^ r {x) 
The interaction of mesons and quarks with the electromagnetic field is described by 

fs 



(45) 



Cr = -eAp(x)p{x)-fQtl> r (x) - eA,{x) ]T 

Z i=l,2,4,5 



fsij + 



The term Cf 71 is generated by minimal substitution with 



(46) 



(47) 



/3ij + 



V3 



where Q is the quark charge matrix and /jjfc are the totally antisymmetric structure constants 
oiSU(3). 

The set of counterterms, denoted by 5C str and 5C em , is explicitly given by 

5C str = SCf + SCf + ££f , 
with 

SCf = 4> r (x) (Z @-M r - S(r) - -f°V(r)]ip r (x), 



SC 



str 



SC 



str 



i) r {x) SM ip r (x), 



i=l 



i=4 



and 

SC ern = _ e A^x)$ r ( X ) (Z - 1) ^Qxjj r ( X ). 
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Here, Z = diag{Z, Z, Z s } is the diagonal matrix of renormalization constants (Z for u, d- 
quarks and Z s for s-quark). The values of Z and Z s are determined by the charge conser- 
vation condition. The simplest way to fix Z and Z s is on the quark level. The same set of 
values for Z and Z s is also obtained when requiring charge conservation on the baryon level. 
Results for Z and Z s will be discussed below. 

Now we briefly explain the role of each counterterm and why the set of constants Z and 
Z s is identical in b~Cf r and 5C em . The counterterm SC em is introduced to guarantee charge 
conservation. The counterterm 5£{ tr , containing the same renormalization constants Z and 
Z s as in S£ em , is added to fulfil electromagnetic local gauge invariance on the Lagrangian 
level. The same term also leads to conservation of the vector current (baryon number 
conservation). Alternatively, 5C em can also be deduced from 5Cf r by minimal substitution. 
In covariant theories the equality of the renormalization constants in 5C{ tr and 5C em is known 
as the Ward identity. The counterterms 5C s 2 tr and 5Cf r compensate the contributions of 
the meson cloud (Fig. la) and meson exchange diagrams (Fig. lb) to the nucleon mass m r N 
(The contribution of meson cloud and exchange diagrams is already taken into account in 
the renormalized quark Lagrangian CZ,.). 

3. Renormalization of nucleon mass and charge 

Now we illustrate the explicit role of the counterterms when performing the calculation of 
the nucleon mass and the nucleon charge. The renormalized nucleon mass rrf N is defined by 
the expectation value of the Hamiltonian TC^ (as derived from the Lagrangian C^) averaged 
over state \<f> > and projected on the respective nucleon states: 

m r N = "< O | / 5(t)d A xH r ^x) |0o > N , (48) 

By inclusion of the counterterms the strong interaction Lagrangian Cf r should give a zero 
contribution to the shift of the renormalized nucleon mass at one loop, that is 

2 -n . 

Am r N = N < O | E "7 / iS (ti) d 4 xi ■ ■ ■ d A x n T[Cf{x x ) . . . Cf{x n )) |0 O >f (49) 

n=l n - J 

= N <M~ll $(ti) d 4 x, d 4 x 2 T[Cf(x 1 )Cf(x 2 )] |0o >c 

3 

- N < O | / 5{t)d A x Y.5£f{x) |0o > N = 0. 

i=i 

The propagator of the renormalized quark field ip r is given by 

iG v (x,y) =< o |T{^ r (x)^ r (y)}|0o u r (x)u r (y) exp[-i£^(x - y o )}0(x o - y Q ). (50) 

It differs from the unperturbed quark propagator iG^(x,y) by terms of order m r , which in 
turn only contribute to the two-loop calculations. Thus, to the order of accuracy we are 
working in (up to one-loop perturbation theory) it is sufficient to use the unperturbed quark 
propagator iG^(x,y) instead of the renormalized one. 

To prove Eq. (f49D, we first note that the contribution of the counterterm b~Cf r is equal 
to zero due to the equation of motion (p8|) , that is 
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^< O | / S^^xSCfix) |0 O > N = 0. (51) 

The counterterms 8Cf r and 5Cf r compensate the contribution of the meson cloud (Fig. la) 
and exchange diagrams (Fig. lb), respectively, with 



- l -J 5(h) d* Xl d A x 2 T[£f(x 1 )Cf r (x 2 )} |0 O >f (52) 
- N < O | / 6{t) d 4 x [5Cf(x) + 5Cf(x)} |0 O > N = 0, 



hence Eq. ( |4"D| ) is fulfilled. The calculation of the nucleon mass m r N at one-loop can then ei- 
ther be done with the " unrenormalized" Lagrangian £ e // or with the "renormalized" version 
£f u ii (|37|) . Both results for rrf N are identical and are given by Eq. (|23|) . 



Now we consider the nucleon charge and prove that the properly introduced counterterms 
guarantee charge conservation. Using Noether's theorem we first derive from the renormal- 
ized Lagrangian fl3~T| ) the electromagnetic, renormalized current operator: 

j? = j$r+j£ +6j» r . (53) 

It contains the quark component , the charged meson component j£ and the contribution 
of the counterterm Sj^r : 

j» r = = \[2 u r ^u r - d r ^d r - s r ^s r } } (54) 



3 

and 

1 

3 



8j$ r = i) r (Z - 1) ^Qi) T = -[2(Z-l) u r ^u r -{Z-l) d r ^d r - (Z s - l)s r ^s T 



The renormalized nucleon charge Q r N at one loop is defined as 

2 -n . 

Qn = 4>ol £ - / S{t)d 4 X d 4 x 1 ...d 4 x n T[Cf(x 1 )...Cf{ Xn )j r {x)} |0o >c ■ (55) 

n=0 n - J 

Charge conservation requires that the nucleon charge is not changed after renormalization, 
that is 

_ J 1 for N = p (proton) . . 

Vn-Vn-^ for N = n (neutron) . W 

Thereby, Qn is the nucleon charge in the three-quark core approximation, which is defined 
as the expectation value of the quark charge operator = J d 3 x j^(x) taken between the 
unperturbed 3q-states |0o >: 



Qn = N < 0o| J 5(t)d 4 xjl(x) |0o >^ . (57) 
Eqs. (p5|) -(p^) completely define the charge conservation condition within our approach. 
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From nucleoli charge conservation at one loop we obtain a condition on the renormal- 
ization constant Z . To fix the constant Z s we should consider the charge conservation of 
baryons containing s-quarks, e.g. S + - baryon. In the one-loop approximation following 
diagrams contribute to the nucleon charge Q r N (see Figs.2a-1): three-quark diagram (Fig.2a) 
with insertion of the quark current j!t r , three-quark diagram (Fig.2b) with the countert- 
erm 5j^ r (three-quark counterterm diagram), meson-cloud diagram (Fig.2c) with the meson 
current j£, vertex correction diagram (Fig.2d) with the quark current j^ r , self-energy di- 
agrams (Figs.2e and 2f) and exchange current diagrams (Figs.2i and 2j) with insertion of 
the quark current j^ r . We also obtain a set of diagrams (Figs.2g, 2h, 2k and 21) generated 
by the counterterms SC^ix) and 5Cf r (x). The contribution of the counterterm S£f r (x) is 
equal to zero due to the equation of motion fl28l) . By definition of the counterterms <5£|* r (x) 
and 5Cl tr (x), the self-energy and the meson exchange current diagrams of Fig.2e, f, i, j are 
compensated by the counterterm diagrams of Fig.2g, h, k and 1, respectively. 

The contribution of the three-quark diagram (Fig. 2a) to the nucleon charge is trivially 
given by 



= N < O | J 5(t)d*xf r (x) |0o > N = Qn- (58) 
The three-quark counterterm diagram (Fig. 2b) is simply related to the one of Fig. 2a with: 
Q r * = {Z-1)QJ = (Z-1) N <4> \ J 5(t)d A xj° r (x)\<f> > N =(Z-l)Q N . (59) 

The meson cloud diagram (Fig. 2c) generates the term 

97 z „ v. 2 °? „*; c 



where 



^ for N = p f ^ for N = p 

n:c I K\c j 

Qn = \ , Qn = < 

-§ for N = n [ | for N = n 

The contribution of the vertex correction diagram (Fig. 2d) is given by 



/^^j;.!?! (61) 



&=iT,K,ri 



where 



i f 2 

for N = p — tj for N = p 




n;d I K-.d I md 

Qn = \ Qn = < Qn = 

| for N = n { -| for N — n 

To guarantee charge conservation, the sum of meson-cloud and vertex correction diagrams 
must be compensated by the counterterm contribution: 
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Q$ + Q r j? + Q r N=U- (62) 

The last requirement fixes the value of the renormalization constant Z at one loop to 

* = 1 - So&fJ^ F ^ (p2) {<py + + ^m}- (63) 

In the two-flavor picture, that is when we restrict to the pion cloud contribution only, we 
obtain a value of Z = 0.9 ± 0.02 for our set of parameters. The contribution of kaon and 
?7-meson loops to the constant Z is strongly suppressed due to the energy denominators in 
Eq. (^3]). In the three-flavor picture we get Z = 0.88 ± 0.03, which deviates only slightly 
from the two- flavor result. The minor role of kaon and 77- meson loop contributions to nucleon 



properties was also found in our previous analysis of meson-nucleon sigma-terms |H|. As 
already mentioned, the renormalization constant Z s is fixed from the charge conservation of 
baryons containing strange quarks (e.g. S + -baryon). Here we obtain the analytical result: 

2 00 

z - = 1 - w &) F ^{^m + 5=M- < 64 » 

In the SU(3) flavor symmetry limit (m u = ma = m s ) both renormalization constants Z and 
Z s are degenerate. Again, charge conservation within our approach is fulfilled both on the 
quark level (when we directly calculate the charge of u, d or s-quark at one loop) and on 
the baryon level. With the value of Z being close to unity for our set of parameters the 
perturbative treatment of the meson cloud is also justified. 



III. ELECTROMAGNETIC FORM FACTORS OF THE NUCLEON 



It is well-known that noncovariant hadron models should explicitly address the question 
of local gauge invariance of electromagnetic form factors (for a detailed discussion see Refs. 
EH" [HI)- Particularly, even when starting from a gauge-invariant Lagrangian the resulting 
physical amplitudes are not necessarily gauge invariant in an arbitrary frame due to the use 
of noncovariant techniques (e.g. noncovariant particle propagators). The lack of Lorentz 
covariance in our approach is linked to the static potential solutions for the quark fields and, 
as consequence, to the use of a noncovariant quark propagator. However, when restricting 
our kinematics to a specific frame, that is the Breit frame, gauge invariance is fulfilled due to 
the decoupling of the time and vector components of the electromagnetic current operator. 

Next we follow the original work by Sachs et al. p4l , |25| and Thomas et al. [^26|27 to set 
up the formalism for the electromagnetic form factors of the nucleon. In the Breit frame the 
initial momentum of the nucleon is p = (E, —q/2), the final momentum is p' = (E, q/2) and 
the 4-momentum of the photon is q = (0, q) with p' = p + q. With the space-like momentum 
transfer squared given as Q 2 = —q 2 = q 2 , in the Breit frame the nucleon charge and 



magnetic (Sachs) form factors are defined by |27 



\J°(0)\N s 



>=X ] n„,Xn s G n e {Q 



(65) 
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Here, J°(0) and J(0) are the time and space component of the electromagnetic current 
operator (|53|); xn s and xjv , are the nucleon spin w.f. in the initial and final state; <jn is the 
nucleon spin matrix. Electromagnetic gauge invariance on the baryon level is fulfilled in the 
Breit frame. In particular, for on-mass shell nucleons we verify the Ward-Takahashi identity 
without referring to the explicit form of the nucleon charge and magnetic form factors: 



+<r-<N4±)\j r (0)\N 8 



■?l> 



-? > 



=0 

Xn. 1 ■ 



<N sl (±)\f r (0)\N s 



■?l> 



(66) 



%a N x q 
2m n 



XNs g n m {Q 2 ) = o 



At zero recoil (q 2 = 0) the Sachs form factors satisfy the following normalization conditions 
G P E (0) = 1, G™(0) = 0, G^(0)=/ip, G n M (0)=^ n , (67) 



where \i v and \i n are the magnetic moments of the proton and the neutron, respectively. 
The charge radii of the nucleons are given by 



< r >\ 



-6 



dG p E (Q 2 ) 



dQ* 



Q 2 = 



< r A > n E = -6 



dG n E {Q 2 



dQ 2 



Q 2 = 



(68) 



< r" > 



N 



6 dG%(Q 



GZ(0) dQ 2 



Q 2 = 



In the PCQM the charge and magnetic form factors of nucleon are defined by 



2 i n f 

Xn ,Xn s G e (Q 2 ) = N < 0o I ~1 / S(t)d 4 xd 4 x 1 . . . d 4 x n e~ 

n=o nl J 



iqx 



(69) 



t ia N x q 



xT[Cf( Xl ) ...£f (x n )j r u (x)]|0 o > 

2 f „ 



AT 
c ' 



2m 



X^(Q 2 ) = O | E \ / * (*) • ■ • d ' x n e~ iqx 

n=0 n - ^ 



xT[£f(x 1 )...£f(x n )>(x)]|0 o > 



IV 



The diagrams of Figs.2a-d contribute to the charge form factor of the nucleon. For the 
magnetic form factors we have to include an additional term due to the "meson-in-flight" 
diagram indicated in Fig. 3. The other possible diagrams at one-loop (like self-energy current 
and meson-exchange current diagrams of Fig. 2) are compensated by the counterterms (see 
previous discussion about renormalization of the nucleon charge). 

In the following we indicate the analytical expressions for the relevant diagrams: 

a) Tree-quark diagram (3g) (Fig. 2a): 



G p E (Q 2 ) 
G p M (Q 2 ) 



3q 
3<z 



G p E (Q 2 ) 

g p m(Q 2 , 



LO 

3q 
LO 

3q 



+ G p E (Q 2 ) 
+ G p M (Q 2 



NLO 

3q 
NLO 

j 

3q 



G n E {Q 2 ) 
GUQ 2 ) 



= 0, 



3q 
3q 



(70) 



-G p M {Q 2 ) 



3q 
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where G N E M (Q 2 ) 



LO 



3q 



are the leading-order (LO) terms of the tree-quark diagram evaluated 



NLO 



is a correction due to the modification 



3q 



with the unperturbed quark w.f. Uo(x); G P E (Q 2 ) 

of the quark w.f. u (x) — ► u^x^rrf) referred to as next-to-leading order (NLO): 

LO 



G P E (Q< 
G P E (Q 2 ) 



3q 
NLO 



Q' 2 " 2 )( l p 2 Q 2 R 2 



3q 



G p M (Q 2 ) 



LO 



4 

Q 2 R 2 
exp[ — | ///' 

12 D 2 



1 + |p 2 4 

Q 2 R 3 p fl + 7p 2 + fp 4 Q 2 R 2 



(71) 



4(1 + \p 2 ) 2 



1 + |P 2 



cxp 



G P M (Q 2 ) 



3q 
NLO 

3q 



Q 2 R 2 \2m N pR 



4 yi + fp 2 ' 



G p M (Q 2 ) 



LO 



m 



3q 1 + \P 2 



Rp (Q 2 R 2 2 - \p 2 



1 + tP 2 



b) Three-quark counterterm (CT) (Fig.2b): 



G P M 



CT 



(Z - 1) G P E {Q 2 



G P m{Q 2 ) = (Z-l)G p M (Q 2 ) 

CT 



LO 

i 

3q 
LO 

3q 



G n E (Q 2 ) 



CT 



o, 



(72) 



) Gm(Q 2 ) 



= ~oG p M (Q^ 

CT O 



CT 



c) Meson-cloud diagram (MC) (Fig. 2c): 

G N E (Q 2 



MC 



G M (Q 2 ) 
where 



4^0 / d *(P 2 +pjQ 2 x)^NN{p 2 ,Q 2 ,x)tUp 2 ,Q 2 ,X 

-1 

oo 1 

y dpp A J dx(l-x 2 ) J r 7T NN(p 2 ,Q 2 ,x)t M (p 2 ,Q 2 ,x) 



3 (g A \ 
MC -m mN {^f) 



MC 



(73) 



F«nn(p 2 , Q%x) = F nNN (p 2 ) F„ NN (p 2 + Q 2 + 2pjQ 2 x), 



t p E (p 2 ,Q 2 ,x) 
t n E (p 2 ,Q 2 ,x) 
t p M (p 2 ,Q 2 ,x) 

t n M (p 2 ,Q 2 ^) 

D^(p 2 ,Q 2 ,x) 

c^ n2 (p 2 ,Q 2 ,x) 



MC 



MC 



MC 



MC 



Cl 1 (p 2 ,Q 2 ,x) + 2C 1 K 1 (p 2 ,Q 2 ,x), 

-Cl 1 (p 2 ,Q 2 ,x)+C 1 K \p 2 ,Q 2 ,x), 

Dl 2 (p 2 ,Q 2 ,x) + *D 22 ( P 2 ,Q 2 ,x), 
5 

-d?(p 2 ,q 2 ,x)-\d 22 (j> 2 ,q 2 ,x), 



2^V,Q 2 ^) 

w$(p 2 ) + w*(p 2 + Q 2 + 2pVCFx) ' 
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d) Vertex-correction diagram (VC) (Fig. 2d): 

LO 

zq ' 200 



Ge(m){Q 2 ) 



vc 



G P E(M)(Q 2 , 



O 2 00 

So (&) I dpp ' ^™(p 2 )^m»(p 2 )|, . < 74 > 



where 





vc 




vc 




vc 


t n M (p 2 ) 


vc 


W*{p 2 ) 



l W «{p 2 )-W K {p 2 ) + ^W v {p 2 ), 

W n (p 2 )-W K (p 2 ), 

^(p 2 ) + l -W K {p 2 ) - ^W V ( P 2 ), 

~l W ^(P 2 ) + 1 W k(p 2 ) + ^W v (p 2 ), 

1 



w%{p 2 )' 

e) Meson-in-flight diagram (MF) (Fig. 3): 



G p E (Q 2 ) 

G P M (Q 2 ) 



= 0, G n E (Q 2 ' 



MF 



MF 



MF 



= 0, G n M (Q 2 ) 



MF 



= ~G P M (Q 2 



MF 



9 

Too 



(75) 



m N 



2 00 1 



We do not include the modification of the quark w.f. in the calculation of one-loop dia- 
grams (Figs. 2b, 2c, 2d and 3) since it only gives corrections to the next (two-loop) order in 
perturbation theory. 

We start our analysis with the magnetic moments of nucleons, ji p and /i n , given by the 
expressions (in units of nucleon magnetons) 



1 + 5- 



+ 



m N / g A V 



LO 



50 \ivFj 
2 

-3"' 

% 
50 



400 WW 



dpp 4 ^atjv(p : 



11 

+ 



26 16 4 



vj3(p 2 ) Wk(p 2 ) U!^(p 2 ) 



K(p 2 ) w 4 k{p 2 ) 

2 00 

1 + *-io(§) 
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wl(p 2 ) w'xip 2 ) w^(p 2 ) 



2 00 



where 



/ip L ° = G p M (0) 



LO _ 2m N pR 
3q ~ 1 + |p 2 



(76) 
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is the leading-order contribution to the proton magnetic moment. The factor 



5 = —rh r R p ■ 



Ip 2 



1 + f P 2 ) 2 



(77) 



defines the NLO correction to the nucleon magnetic moments due to the modification of 
the quark w.f. (see Eq. ([71])). Note that the well-known SUq relation between nucleon 
magnetic moments fi n /fi p = —2/3 can be easily deduced from Eq. ( [7E| ) if: i) we restrict 
to contributions from one-body diagrams in Figs.2a-d, corresponding to the additive quark 



picture, and ii) apply the S77(3)-fiavor limit (Mo- 
have 



M, 



K 



M„ 



(SUe) = _±.XSU 6 ) _ LO 
H>p — nPn — H-p 



6 M 



+ 50 mN 



9a 



dpp 



F 2 

4 nNN 



200 

(p 2 



9a 

7lF 



dpp 



M$). In particular, we 
(p 2 ) 



rp2 
4 r ivNN 



wl(p 2 



w%(p 2 



(78) 



Taking into account the meson- in-flight diagram (Fig. 3) generated by two-body forces leads 
to a deviation of the ratio \x n j \i v from the naive SU 6 result. 
In general case, for our set of parameters we obtain: 



p p = 2.62 ± 0.02, fj, n = -2.02 ± 0.02, and 



P P 



-0.76 ±0.01 



(79) 



where the range of theoretical predictions reflects the variation of the size parameter R. The 
separate LO (Pp°, Pn°) an d NLO (fip LO , yU^ i0 ) contributions to the magnetic moments 
are given by 



1.8 ±0.15, 
= ji p — ji p 



LO 



-^P 



LO 



(80) 



0.82 ±0.13, 



NLO 



Pn - yU 



LO 



-0.82 ±0.08. 



Our results for \x v and //„ are close to the ones of the chiral soliton quark model by Diakonov 
and Petrov 0: \i p = 2.98 and Li n = —2.26. For the electromagnetic nucleon radii we obtain 



M 



0.84 ± 0.05 fm, 
: 0.82 ± 0.02 fm, 



< r 2 >»= -0.046 ± 0.006 fm 2 , 



(81) 



M 



0.85 ± 0.01 fm. 



The LO contributions to the charge radius of the proton (see Eq. (0)) and to the magnetic 
radii of proton and neutron are dominant 



r p 
r E 



LO 



0.77 ± 0.06 fm, 



M 



= r 



LO 



M 



0.73 ± 0.06 fm. 



LO 



(82) 



For the neutron charge radius squared we get the observed (negative) sign, but its magni- 
tude is smaller than the experimental value. As in the naive SU(6) quark model, the LO 
contribution to the neutron charge radius is zero and only one-loop diagrams give nontrivial 
contributions to this quantity: the meson-cloud (MC) < r 2 >e-mc = ~ 0.072 ± 0.006 fm 2 
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and vertex-correction (VC) < r 2 >£.yc~ 0.026 fm 2 diagrams. The total contribution to the 
neutron charge radius is given by 

< r 2 > n E =< r 2 >% MC + <r 2 > n E . vc = -0.046 ± 0.006 fm 2 . (83) 

The magnitude of our result for the neutron charge radius, which is too small compared 
to the experimental value, is due to the reduced contribution of the pion cloud diagram 
(Fig. 2c). The size of the meson cloud effects is in turn related to the magnitude of the 
renormalization constant, which for our model is Z=0.9 and therefore close to unity. With an 
increased contribution of the pion cloud, that is a decrease of Z, one is able to reproduce the 
experimental result for the neutron charge radius, but one also obtains a worse description 
of the magnetic moments and the proton charge radius. Also, perturbation theory up to one 
loop becomes less reliable in terms of good convergence properties. For example, choosing 
p = 0.12 (or g A = 1.5) and R = 0.52 fm we obtain < r 2 > n E = -0.08 fm 2 and Z = 0.73. 
However, for the other quantities we get a worse fit with: /i p = 2.56, fi n = —2.18, r p E = 0.75 
fm, r p M = 0.9 fm and = 1 fm. In comparison, the cloudy bag model gets a value of 
< r 2 > E = -0.14 fm 2 by an increased contribution of the pion cloud with Z m 0.73 ||27|| . 
As a result, the predicted values for the proton magnetic moment are also less satisfactory 
with fi p = 2.2 as in the original version of the cloudy bag model || and fi p = 2.46 in the 
improved version |p7|| . The prediction for the proton charge radius is also relatively small: 



r p E = 0.73 fm §. 

An improved fit to the neutron charge radius, while keeping the good results for the other 
static observables, can possibly be achieved by including center-of-momentum corrections 
and also Lorentz boost effects Refs. P,|5|,|3O|-07[|), which are based on approximate techniques 
and are not included yet in the present evaluation. Furthermore, additional contributions of 
excited quark/antiquark states in the meson loop diagrams will further enhance the meson- 
cloud effects, which up to now were not studied yet consistently. 

In Table I we summarize our results for the static electromagnetic properties of the 
nucleon in comparison with experimental data [ 38| . Finally, in Figs. 4-8 we indicate the Q 2 - 
dependence of the electromagnetic form factors. In Fig.4 we present the result for the proton 
charge form factor at a typical value for the free size parameter R = 0.6 fm. The dotted 
line corresponds to the LO result, the short-dashed line marks the NLO contribution, while 
the solid line corresponds to the total result at one loop. And, finally, the long-dashed line 
indicates the dipole fit to experimental data given by 

Gd{Q2) = (l + QVO^lGeV 2 ) 2 ' (84) 

The dependence of the proton charge form factor on the free parameter R is presented in 
Fig. 5. We plot three curves corresponding to R = 0.55 fm, R = 0.6 fm and R = 0.65 fm. 
An increase of R leads to a decrease of the proton charge form factor and vice versa. In 
Fig. 6 we give the results for the proton charge and the nucleon magnetic form factors for 
R = 0.6 fm and compare them with the common dipole fit. In Fig. 7 we plot our predictions 
for the ratio fi p G E (Q 2 ) / G P M (Q 2 ) at R = 0.6 fm and compare them with experimental data 
taken from |3{| and [flOX Our predictions for the ratio G^iQ 2 ) j \\x n G d{Q 2 )) are presented 



in Fig.8. Experimental data are taken from [^TJ- |44|| . Finally, in Fig. 9 we plot the results 
for the neutron charge form factor at different values of R: R = 0.55 fm, R = 0.6 fm and 



R = 0.65 fm. Our results are lower than the experimental points taken from Refs. |45|- ||49 
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IV. SUMMARY AND CONCLUSIONS 



In this paper, we have considered the perturbative chiral quark model (PCQM) based 
on an effective chiral Lagrangian which includes confinement phenomenologically. The La- 
grangian basically describes nucleons as bound states of three valence quarks surrounded by 
a perturbative cloud of pseudoscalar mesons as dictated by the chiral symmetry requirement. 

The main aim of this investigation was to present a consistent formalism when treating 
the Lagrangian in the one-loop expansion. We thereby employed a technique, which by 
appropriate introduction of counterterms allows consistent mass renormalization both on the 
quark and the nucleon level. In a further extention of the PCQM we considered the coupling 
of the photon field to the nucleon. Again, charge renormalization was introduced to keep 
a proper definition of the nucleon charge. Local gauge invariance of the electromagnetic 
interaction is fulfilled on the Lagrangian level by construction. Due to the noncovariant 
nature of the effective confinement we introduced, local gauge invariance is not necessarily 
fulfilled for physical amplitudes in any reference frame. Only when working in the Breit 
frame, nucleon matrix elements were shown to be consistent with the Ward identity without 
further need to introduce extra terms, like contact interactions, in the Lagrangian. 

As one basic application of the developed formalism we considered the electromagnetic 
form factors of the nucleon. Two approximations were introduced: 

i) restriction of the relativistic quark states to the ground state, that is N and A inter- 
mediate states occur in the one-loop terms; 

ii) Gaussian form of the ground state single quark wave function, modelled for low-Q 2 
physics of the calculated observables. 

The derived quantities (magnetic moments, radii and form factors) contain only one 
model parameter R, which is related to the radius of the three-quark core, and are otherwise 
expressed in terms of fundamental parameters of low-energy hadron physics: axial coupling 
constant gA, weak pion decay constant F and set of QCD parameters (current quark masses 
rh and m s and quark condensate parameter B). Predictions are given for a variation of the 
free parameter R in a physical region from 0.55 fm to 0.65 fm corresponding to < r\ >f^_ core 
ranging from 0.5 fm 2 to 0.7 fm 2 . Our low Q 2 results, that is for the nucleon magnetic 
moments, the charge radius of the proton and the magnetic radii of the nucleons are in 
reasonable agreement with data, strengthening the phenomenological validity of the PCQM. 
For the charge radius of the neutron we obtain the correct sign but its magnitude is smaller 
than the experimental value. 

We also presented a detailed discussion of the Q 2 -dependence of the nucleon form factors 
in the space-like region. Although the model underestimates the finite Q 2 behaviour of 
the form factors, certain additional physics aspects, such as cm. correction and Lorentz 
boost effects, are known to improve the phenomenological fit (see detailed discussion in 
Refs. P,|5|,|5D|-p7|] ) . The issue of cm. corrections as evaluated in the PCQM machinery 
and their consistency with the chiral expansion of the nucleon mass will be addressed in 
a forthcoming publication. With the formal and phenomenological foundation established, 
the proposed model serves as a basis for further applications to recent issues of nucleon 
structure, for example such as photo- and electroproduction of pseudoscalar mesons, which 
will be pursued in future. 
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TABLES 



Table I. Static nucleon properties. 



Quantity 


Our Approach 


Experiment 38 


Hp 


2.62 ± 0.02 


2.793 




-2.02 ± 0.02 


-1.913 


ym 


-0.76 ± 0.01 


-0.68 


r E (fm) 


0.84 ± 0.05 


0.86 ± 0.01 


< r 2 > E (fm 2 ) 


-0.046 ± 0.006 


-0.119 ± 0.004 


r p M (fm) 


0.82 ± 0.02 


0.86 ± 0.06 


r n M (fm) 


0.85 ± 0.01 


0.88 ± 0.07 
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FIGURES 



Fig.l: Diagrams contributing to the baryon energy shift: 
meson cloud (la) and exchange diagram (lb). 



Fig. 2: Diagrams contributing to the nucleon charge: 

triangle diagram (2a), triangle counterterm diagram (2b), meson-cloud diagram (2c), ver- 
tex correction diagram (2d), self-energy diagrams (2e) and (2f), self-energy counterterm 
diagrams (2g) and (2h), exchange current diagrams (2i) and (2j) and exchange current 
counterterm diagrams (2k) and (21). 



Fig. 3: Meson-in-flight diagram. 

Fig. 4: Proton charge form factor: LO, NLO and Total contributions. 



Fig. 5: Proton charge form factor for different values of the 3q-core radius R = 0.55, 0.6 and 
0.65 fm. 



Fig. 6: Proton charge and nucleon magnetic form factors at R = 0.6 fm. 



Fig. 7: The ratio fi p G p E (Q 2 ) /G P M (Q 2 ) at R = 0.6 fm. Experimental data are taken from |5S 
(JLAB) and Q (MIT-Bates). 



Fig. 8: The ratio G^^Q 2 ) j \\i n G d{Q 2 )) at i? = 0.6 fm. Experimental data are taken from 
|H (MIT-Bates), g| (NIKHEF), || (ELSA) and @ (MAMI). 



Fig. 9: Neutron charge form factor for different values of R = 0.55, 0.6 and 0.65 fm. Experi- 
mental data are taken from g§| (MAMI-1), g§| (MAMI-2), (MAMI-3), g§] (MIT) and 
(PARIS). 
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